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■ Abstract 
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! It is pointed out that the retardation terms given in the original 

I Fermi-Breit potential vanish in the center of mass frame. The re- 

■ tarded one-gluon exchange potential is rederived in this paper from the 
"^H , three-dimensional one-gluon exchange kernel which appears in the ex- 

CL|| act three-dimensional relativistic equation for quark-antiquark bound 

Q^' states. The retardation part of the potential given in the approxima- 

D . tion of order /m^ is shown to be different from those derived in the 

previous literature. This part is off-shell and does no longer vanish in 

' the center of mass frame, 

[x^ '. PACS number(s): 12.40 Qq, ll.lO.St, 14.40.Gx. 

^ \ Key words: One-gluon exchamge potential, Retardation effect. 

In the quark potential model'^"^', the one-gluon exchange potential (OGEP), 
or say, the Fermi-Breit (F-B) potential, usually is chosen to describe the 
short-range interaction between two quarks (gg), two antiquarks (gg) or a 
quark and an antiquark {qq). This potential was originally derived from the 
t-channel one-gluon exchange S-matrixt''' or the t-channel one-gluon exchange 
kernel appearing in the Bethe-Salpeter (B-S) equationt*^'^! . In the relativistic 
version, the potential may be represented in the form 

V{m,;PuP2) = U+{Y,)U+{Y^)V{l,2)U{-[^,)U{-[^2) (1) 
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where 

is the Dirac spinor for ith quark, satisfying the normahzation condition 
U+{^i) U{^i) = 1 and 



V{1,2) = 27i72^(l,2) (3) 



is the interaction operator. The first equahty in Eq.(3) shows the general 
relation of the interaction operator with the B-S irreducible kernel K(l,2), 
while the second equality gives the expression for the t-channel one-gluon 
exchange interaction in which Cg stands for the color matrix 



-T^-, for qq. 

and D^iy{Q) denotes the gluon propagator whose argument Q is defined as 

<9 = p'l - Pi = P2 - p'2 (5) 

Upon inserting Eq.(3) into Eq.(l), in the Feynman gauge which was often 
used in the previous derivation of the FBP, we may write 

V{p'i,p'2;PuP2) = ^nasCsTCp"'^, "p'2; P^i, 'P2)D{Q) (6) 
where as — 5'^/47r, 

nVi, 'p'2; 'pi: 'P2) = u^{VMriU{-^i)u^{VMi2,u{-^2) (7) 

and 

Ordinarily, the OGEP is derived in the nonrelativistic approximation up to 
the order of j?- j nx^ [v^ j (?) and given in the three-dimensional space. By 
making use of the nonrelativistic expression for the quark energy 



the three-dimensional spinor matrix in Eq.(7) will be written, up to the order 
/m? , as foUows'^'^l 



Smyrna 

mi — mo -r^s 0^ A; 



2mim2mi2 mim2 

+- [P X Y - A; ■ — ^ 

4mi2 mi m2 

^ ^ ' ^1-^2 + -Ti2{^ - k ) 



6mim2 Amim2 

+T^(^ X ^) • [(2 + — + (2 + — )^2] 
4mim2 mi m2 



where 



7i2(T - A; ) = ^ "'^ ^ 2 "^! - k )~a2 ■ {"q - k) - ^~ai ■ ~a2 (H) 
~q — k 

is the tensor force. In the above, we have introduced the total momentum 
P and the relative momentum k and ~q for the initial and final states of 
the two particles respectively. They are defined by 

p^i = rjiP + k ,J)'2 = 1I2P - k (12) 

P^'i = r]iP +~q,~p'2 = r]2P -~q 

rrii . -, r. 

r)i = ,mi2 = mi +m2,^ = 1,2 

mi2 

To obtain a three-dimensional potential, the question left is centered on how 
to reduce the four-dimensional gluon propagator to a three-dimensional form. 
In the original derivation of the F-B potential, with the recognition that 
the Coulomb term in the potential is dominant 1^'^^, meaning that the ratio 
Qq / Q ^ is much less than unity, the vector boson propagator was reasonably 
expanded as'^^ 

D{Q) = ^ ^ + #) (13) 
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Based on the energy conservation denoted in Eq.(5) and the nonrelativistic 
fermion energies written in Eq.(9), in the approximation of the order /m? 
, the time-component of the transfer momentum squared in Eq.(13) may be 
taken to be 

Ql = iE[ - E,)iE2 - E',) (14) 

~ -^—{uj ■ ^1 g ■ ^2 + 2g ' g ■ ^2 - 2 g ' g • - g ^) 

4mim2 

On substituting this expression into Eq.(13), the propagator was approxi- 
mately represented in the following three-dimensional form^^'^l 

D(g) = ^-^— + (15) 

Q2 4mim2 2mim2Q^ mxm^Q^ 

where the first term gives the coulombic potential which describes the in- 
stantaneous interaction, while the remaining terms represent the retardation 
effect. When Eqs.(lO) and (15) are inserted into Eq.(6) and the terms higher 
than p^/ VP? are neglected, the well-known F-B potential will be recovered. 

It should be pointed out that according to the procedure of derivation 
stated above, the retardation terms in the F-B potential may possibly be 
given different expressions. This is because the Q\ In Eq.(13) may be repre- 
sented in different ways. For example, in accordance with the relation shown 
in Eq.(5), we may also write Q\ = {E[ — Ei)^ or Ql = {E2 — -^2)^. Cor- 
responding to either of these two representations, in the approximation of 
order p^/m^, the propagator will be given a different expression 

g2 4m? ^2g4 ^2Q2 

or 

g2 4mi rnlQ^ mlQ^ 

The different expressions given in Eqs.(15)-(17) will be consistent with each 
other when and only when they are constrained by the energy-momentum 
conservation. In fact, when the nonrelativistic energies in Eq.(9) are substi- 
tuted into the energy-conservation 

E[ + E'^ = Ei + E2 (18) 
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and then the momentum conservation shown in Eq. (5) is employed, one may 
find 

— >• 2 

= {rni~p 2 - rn2~p i) ■ Q (19) 

With this relation, it is easy to prove that the expressions in Eqs.(15)-(17) 
may be transformed to one another. This proof indicates that the energy- 
momentum conservation acts a necessary constraint to guarantee the retar- 
dation terms in the F-B potential to be uniquely determined. 

When the energies in Eq.(18) are replaced by the nonrelativistic ones 
shown in Eq.(9) and then the momenta of single particles are represented 
through the total and relative momenta defined in Eq.(12) , it will be found 

= (20) 

With the aid of this relation and those shown in Eq.(12), the Ql defined in 
the first line of Eq.(14) can be rewritten as 

Ql = ^y^-ri)iri-V2^) (21) 



This result may directly be derived from the last expression of Eq.(14), when 
the momenta 'pti and in the expression are replaced by the total and 
relative momenta and then using the following relation 

Q-~k^ -^Q' (22) 

which is obtained from Eq.(19) by employing the relations denoted in Eq.(12). 
From Eq.(21) , it is clear to see that in the center of mass frame ( P = ), 
we have 

Ql = (23) 
Due to the above equality, Eq.(13) becomes 

D{Q) = X (24) 

which is completely instantaneous. This result shows that the ordinary F-B 
potential fails to describe the retardation effect existing in the center of mass 
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frame. The failure obviously originates from the use of the on-mass-shell 
approximation shown in Eq.(14) in the derivation of the three-dimensional 
approximate propagator given in Eq.(15). The above result also indicates 
that only the instantaneous part of the F-B potential can be reasonably used 
in the study of hadron spectroscopy. 

Quite recently, the retardation effect has attracted some interest in the 
study of hadron spectroscopy^^'^'^l . In Ref.(8), the authors derived a certain 
retardation terms for the scalar and vector couplings in a semirelativistic 
quark potential model. In Ref.(9) , the retardation effect was especially 
investigated in some detail. In the investigation, the authors did not give 
an explicit three-dimensional expression of the retardation corrections to the 
potential. Instead, they directly evaluated the retardation corrections to the 
quarkonium with the help of a regularization procedure. The retardation 
corrections were defined as the difference between the values obtained from 
solving the four dimensional B-S equation for bound states by using the exact 
one gluon exchange B-S kernel and the instantaneous one. 

In this paper, we study the retardation effect from a different angle and 
are devoted to deriving explicit expressions of the retardation part of the 
OGEP. The essential feature of the derivation is that we start from an exact 
three-dimensional relativistic equation for qq{qq or qq ) bound states which 
is equivalent to the corresponding four-dimensional B-S equation'^^'^^l. 



where E is the total energy of the system, Xpa{~q^) is the wave function for 
the system. 



is the Hamiltonian for the ith fermion and K{~p , ~q , k , E) denotes the three- 
dimensional interaction kernel which may be explicitly expressed in terms of 
a few types of Green's functions and is energy- dependent as any interaction 
kernel'^^'^^l. In this paper, we only concern the t-channel one-gluon exchange 
kernel. This kernel still has the form as represented in Eq.(3) except that 
the gluon propagator becomes three-dimensional. In the Feynman gauge, the 




[E - hiCpi) - h{'p2)\Xpcc{'q) 
——^K{'p', ~q, k , E)xpa{ k ) 



(25) 



K'P^i) = "a i • Pi' + Ami, i = 1, 2 



(26) 
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propagator is of the forml^^'^^l 

D(7,r,t,E) = -1^1— ^^--1^^^-^^ (27) 

1 

+ 



E - ECp'i) - ECp"^) - Q + 

where Q = 1^ — k is the magnitude of the three-dimensional transfer 
momentum, E is the total energy of the interacting system, £"(^1), £"(^2), 
E[ p ^) and E[ p 2) are the energies of single free particles which are on the 
mass-shell. It is easy to verify that when the total energy in Eqs.(8) and (27) 
is put on the mass shell, E = E{'J)^i)+ E{'J^2) = E{'p'^) + E{'J^2)j the both 
expressions in Eqs.(8) and (27) may be transformed to each other. When the 
energy E is off-shell, the propagator shown above just is the exact three- 
dimensional representation of the propagator written in Eq.(8) and therefore 
containing all the retardation effects in the case of bound states. 

Now let us discuss the retardation terms in the OGEP. It is clear that 
when the propagator in Eq.(27) is inserted into Eq.(6), we obtain a rigor- 
ous relativistic representation of the retarded OGEP written in the three- 
dimensional space. In the following, we are interested in deriving a retarded 
OGEP in the nonrelativistic approximation which is used in the nonrela- 
tivistic quark potential model. For this purpose, it is necessary to write 
the propagator in Eq.(27) in the nonrelativistic approximation of the order 
p^ /rn? by employing the approximate expression in Eq.(9) 

D(P,^,t) = -L[ (28) 



where e = E — mi — 1712 is the binding energy. Firstly, we discuss a special 
case that the binding energy is small, being of the order of p'^/m^. In this 
case, noticing the common recognition mentioned before for the expansion 
shown in Eq.(13) that in the nonrelativistic approximation, the Coulomb 
term plays a dominant role in the OGEP^'^'^J, the above propagator may 
approximately be taken to be^^^l 
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where fi is the reduced mass of two particles, P , "g*and k are the momenta 
defined in Eq.(12). In Eq.(29), except for the first term which is instanta- 
neous, the other terms all come from the retardation effect. In contrast to 
the corresponding expression given in Eq.(15), the propagator in Eq.(29) is 
off-shell, i.e., energy-dependent, and, as we see, the retardation terms in it 
no longer vanish in the center of mass frame. Substituting Eqs.(lO) and (29) 
into Eq.(6), the OGEP may be separated into two parts: the instantaneous 
part Vi{P ,~q', k) and the retardation part Vr{P ,~q , k) 

V(P, 'k) = Vj(P, 'k) + Vr(p, 'k) (30) 

where 

vi(p, = ^p^r( p , ^) (31) 

— k 

in which the r( P , ~q , k ) was represented in Eq.(lO) and 

47ra,C, , P2 1 
— k 

By the Fourier transformation 

cP'q (fk 

(2^(2^ 
= 5^(1^ -T'')V(P,T') 

We may obtain a local OGEP in the position space 

v(p, r ) = Vi(p, r) + Vr(p, r) (34) 

where V/( P , T^) and Vr( P , T*) are respectively the instantaneous part and 
the retardation part of the potential V{ P , T^). The instantaneous part is'-*^'^^' 

v,(y,r) = (35) 

r mi2 277117712^12?" 

-2r^y . v^) - ^K-^^)V ^3(y) '—Jr^V%. 

277ii77i2 mim2r^ 



y.(P , k) ^ ^^[e - ^ - -(T^ + k ^)] (32) 



V(P,r,r') = /7^7^e^^-^-^^-^V(P,r,^) (33) 
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where 



omim2 4mim2r^ 

-ir^—2^r ■ [(2 + ^)Sl + (2 + 

zmim2r^ mi m2 



P = -^V^, 5=1 = ^,^ = 1,2 (36) 
i-ii = X (-iV-^),^^ = X (-iVy*) 

^(^) = ^"ct 1 • y 0^2 • ^ - ^"^1 • 0^2 



i? and denote the center of mass coordinate and the relative coordinate 
respectively. For the retardation part, we have^^^l 

Vni'P, r) = --asCs{e In - - + -1^[1 + 21^ • + 2r^ In(-) V^^]} 
TT To 2mi2 4:iir^ Tq 

(37) 

where the Fourier transformation 

' mil (38) 



in which ro is a size parameter has been used. The logarithmic function in 
Eq.(37) was ever phenomenologically introduced to the quark potential model 
in some previous works^^^^ for describing the intermediate-range interaction. 
Now, it appears as a natural result of the retardation effect. 

Next, we have a particular interest in deriving a retarded OGEP in a 
more general case that the binding energy is not required to be small. In this 
case, the gluon propagator in Eq.(28) should be expanded in the form in the 
approximation of order /w? as follows 

^^^'^'^^ = g(^^^-g^[£ + i^^^ + ^^)]>^^^) 

= i^ + L»^(P,Y,"fc') 
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where 

P2 1 



Dn{P,r, k) = D,{Q)-D2m7^ + ^i^'+ k')] (40) 

2mi2 4// 



in which 

and 

With the propagator given in Eq.(39), the total OGEP may also be divided 
into the two parts as written in Eq.(30). The instantaneous part of the 
potential is still the one represented in Eq.(31) in the momentum space and 
the one in Eq.(35) in the position space. While, the retardation part of the 
potential now is of the following form in the approximation of order p^/rn^ 

Vr(P, = 47ra,C,{Di(g)r( P , 'kyD^m^^+^it'+'k')]} 

zmi2 

(43) 

Through the Fourier transformation as written in Eq.(33) , we may obtain 
from Eq. (43) a local form of the retardation part of the potential as follows 



rrii 7712 mim2 mim2 

_ lijn^-m^) ^ [Vr A(r) + 2D,{r)V^] (44) 
2mim2mi2 

xV,l.|(2 + ^)S; + (2 + ^)-Kl-|^fl.M 
+ J-[V|.D2(»-) + 2VrD2(r) ■ Vr + 2D^(T)Vy\} 
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where 

Di{r) 



e e 



(45) 



Airr 
1 



cos Xq 



/ 



dx e 



27r2 i x"^ + Xn 
^ 

2 2 
{cos xo — 1 -\ — In xo sin xq (xq) cos xq — vi^o) sin xq] } 



in which xq = \s\ r, 



TT 



TT 



2n(2n)! 



«^»'=S(2.+i)(2.+i)r''<^°) = ^\?, 



(46) 



and 



D2{r) 



(27r) 



/ 



g ( Q -ef 
sinx 



27r2 i (x + xo)2 

1 .TT 



(47) 



27r2^2 



[— sin Xq — In cos Xq — Hxq) sin — rj{xQ) COS Xo] 



In the practical calculation, the integral expressions in Eqs.(45) and (47) are 
more suitable to be used. In the case of the binding energy being small, the 
functions Di[r ) and D2{r) will be approximated to 



Di{r) 



27r2 



ln(|£| r) 



and 



27r2 



[1 + Mkk)] 



(48) 



(49) 



With this approximation, the expression in Eq.(44) will be reduced to the 
one given in Eq.(37) if we set |£| = Tq ^ in the function ln(|£| r). 
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